Quantum chaos and critical behavior on a chip 
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The Dicke model describes N qubits (or two-level atoms) homogenously coupled to a bosonic 
mode. Here we examine an open-system realization of the Dicke model, which contains critical and 
chaotic behaviour. In particular, we extend this model to include an additional open transport qubit 
(TQ) (coupled to the bosonic mode) for passive and active measurements. We illustrate how the 
scaling (in the number of qubits N) of the superradiant phase transition can be observed in both 
current and current-noise measurements through the transport qubit. Using a master equation, we 
also investigate how the phase transition is affected by the back-action from the transport qubit and 
losses in the cavity. In addition, we show that the non-integrable quantum chaotic character of the 
Dicke model is retained in an open-system environment. We propose how all of these effects could 
been seen in a circuit QED system formed from an array of superconducting qubits, or an atom 
chip, coupled to a quantized resonant cavity (e.g., a microwave transmission line). 



I. INTRODUCTION 

Understanding and categorizing complex modes of be- 
havior, such as quantum phase transitions^ and quantum 
chao8^, is an important part of quantum many-body the- 
ory. Recently, concepts and formalisms from quantum 
information theory have been used to understand and 
classify several aspects of criticality^i^i^i^ii. However, the 
realization of strong coupling regimes, coherent dynam- 
ics, and careful readout necessary to observe these phe- 
nomena in laboratory conditions is challenging. 

Our goal here is to show how a particular quantum 
phase transition, the Dicke superradiant transition^^^^^ 
behaves when coupled to the environment and measured 
using transport techniques, as is the case in realistic ex- 
perimental conditions. The Dicke model describes N 
two-level "atoms" or qubits coupled to a common single- 
mode cavity. We focus on this model because of the re- 
cent advances in on-chip 'circuit- QED 'i^^iiii^ii^iiiii^ii^, 
where the strong coupling regime is accessible, and which 
allow for coupling to a range of artificial atoms and mea- 
surement apparatuses. In particular, we propose a dis- 
persive measurement scheme to observe this transition 
by coupling either a superconducting qubit array, or an 
atom chip, to a cavity which is simultaneously (disper- 
sively) coupled to a non-equilibrium measurement device 
(a so-called ''transport" qubit^^), realisable with a super- 
conducting single electron transistor, or double quantum 
dot. The geometry of the proposed device is shown in 
Fig. 1 and described in detail in its caption. 

We begin by outlining the salient features of the phase 
transition in the Dicke model, and existing work in this 
area (section II), and discussing the closed (section HI) 
and open (section IV) descriptions of the model. We then 
investigate how coupling to a transport qubit allows read- 
out of the phase transition properties, and give analytical 
and numerical results for the current and current-noise 



in the zero back-action limit (section V). We then dis- 
cuss our main result: that the current can be used as 
an observable order parameter to detect the phase tran- 
sition (section VI). This complements a recent surge of 
interest in identifying signatures of complex behaviour 
in mesoscopic transport measurements^^^^^^^^. We then 
consider back-action and decoherence (cavity-loss) ef- 
fects using a master equation approach (section VII). 
We show that both transport qubit back-action and cav- 
ity loss appear to only have a weak affect on the cur- 
rent measurement near the critical point. In addition, 
we show that the Liouvillian describing the open-system 
dynamics has an eigenvalue spectrum similar to that of 
the Wigner-Dyson distribution of random matrix theory, 
as in the closed system case (section VIII). Finally, we 
briefly discuss practical schemes to realize this model in 



an experiment' 
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(section IX). 



We point out that the properties we investigate here re- 
quire the precise control of the couplings between qubits 
and the resonator, and access to a very strong coupling 
regime, both of which are difficult to achieve. However it 
was recently shown^l that the generalized Dicke model, a 
variation of the Dicke model where the couplings between 
the N qubits and the cavity are inhomogeneous, still has 
all the critical properties of the standard Dicke model. 
This indicates the universality of the Dicke phase transi- 
tion, as well as making an experimental realization more 
feasible. Furthermore, the critical point in the transition 
relies on the relative values of the coupling strengths and 
the level-splittings of the qubits. Thus, while the qubit- 
boson coupling strengths are not tunable in a real ex- 
periment, the level splittings of the qubits can typically 
be controlled by external parameters (e.g., in the case 
of superconducting flux qubits via an external magnetic 
flux^^). Furthermore, a realization using Raman transi- 
tions in atoms in an optical cavity has been proposed as a 
method to reach a controllable strong coupling regime^. 
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FIG. 1: (Color Online) Geometry of the proposed device. N 
qubits are placed at the anti-node(s) of a 'cavity' with res- 
onant frequency oo (depending on specific realization, they 
might alternatively be placed at one anti-node, or at sub- 
sequent anti- nodes, to minimize qubit-qubit interactions). 
Their energy splittings are on-resonance with the oscillator 
ooo = 00. An additional 'open' transport qubit (TQ, shown in 
the left) is also coupled to the cavity off-resonance (e / a;), 
and used to passively read-out the state of the cavity mode. 
In the figure, the solid lines represent tunneling, while the 
dashed lines represent energy gaps. The properties of the 
transport qubit are defined by an energy splitting e, coherent 
tunneling rate A and transport rates Tl and Tr. In addi- 
tion, the cavity mode has a cavity decay rate 7^, not shown 
in the figure. Superconducting artificial atoms coupled to an 
on-chip cavity (e.g. a quantized LC oscillator or microwave 
transmission line) are a feasible realization using current tech- 
nology. The transport qubit can be realised using a charge 
qubit in the transport regime, i.e., a superconducting single 
electron transistor. Alternatively, a large number of qubits 
in the form of two- level atoms in an atom chip, coupled to 
a transmission line, has recently been proposed as a way to 
realise the large- A/" Dicke modet^. 



II. DICKE SUPERRADIANT PHASE 
TRANSITION. 



Historically, the Dicke Hamiltonian (DH) describes the 
dipole interaction between N "atoms" and bosonic 
field modes. Typically^^ the atoms are considered to be 
at fixed sites within a cavity of volume V. The atoms are 
assumed to be well separated, and thus non-interacting. 
Hereafter we refer to the atoms as 'qubits', and the ad- 
ditional measurement qubit as the 'transport qubit'. To 
observe critical phenomena we consider the single-mode 
case with = 1. We do not make the rotating wave ap- 
proximation, allowing the model to describe both weak 
and strong coupling regimes (and we omit the term). 

Previous work^^^ on this model has shown an exact 
analytical solution in the limit TV ^ 00. Furthermore, 
the transition was characterized as a breaking of parity 
symmetry at a particular value of the coupling between 
qubits and cavity (denoted by A, with the critical value 
being Ac). Both the qubits and the cavity bosonic de- 
grees of freedom become 'macroscopically occupied' (i.e., 
of 0{N), the number of qubits) in the regime above the 
critical point A > Ac- For finite arrays of qubits. A/", the 
system is known to exhibit power-law scaling^, quantum 



chaos^, and critical entanglement 

Several proposals for an experimental realisation of this 
system have already been made. For example, Dimer et 
al^^ proposed a cavity QED realisation, and discussed in 
detail the effect of the cavity decay on the phase transi- 
tion. In another work, Chen et al^^ proposed using su- 
perconducting charge qubits coupled to an optical cavity, 
so that the critical properties can be observed in the op- 
tical mode using heterodyne detection. In addition they 
proposed observing the phase transition as a function of 
level splitting, as discussed earlier. 



III. DICKE HAMILTONIAN. 



The single-mode Dicke Hamiltonian is defined as 
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= uqJz uja^a-^ —={a^ ^ a){J^ ^ (1) 
V N 



where = YliLi^l^ J± = Ylf=i^^± collective an- 
gular momentum operators for a pseudo-spin of length 
j = N/2. These operators obey the usual angular mo- 
mentum commutation relations, [Jz,J±] = ^J± and 
[J+, J_] = 2Jz. The frequency coq describes the qubit 
level splitting, uj is the oscillator field frequency, and A 
the qubit-field coupling strength. Because of their mu- 
tual interaction with the oscillator field the qubits are not 
independent. The A/a/TV scaling is important to realise 
the thermodynamic limit. It essentially bosonises the 
low-energy part of the state space of the collective angu- 
lar momentum. Physically, this scaling implies that the 
density of qubits is constant, so that the cavity volume 
becomes larger as N is increased, consequently reducing 
the electric field density and thus the effective interaction 
with each individual qubit. 

First, we will show analytical results for an en- 
tirely passive measurement of the system, using an off- 
resonance ancillary qubit with current transport. Sec- 
ondly, we will treat the back-action of the ancillary qubit 
as a fully quantum interaction, with Markovian transport 
properties, and including decay terms for the cavity. We 
will see how this alters the final current measurements, 
as well as how it changes the properties of the phase 
transition. 



IV. MASTER EQUATION. 

To take into account the back-action of the transport 
qubit on the Dicke Hamiltonian, we can model the whole 
system using a master equation. 
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Pit) = L[p{t)]=-z[H,p{t)]^Lo[p{t)] 



(2) 



H = Hn + H7 



Lf) = Ln 
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where 



PASSIVE MEASUREMENT. 



is the Hamiltonian of the transport qubit, 



^int = gcr^a'^a 



(3) 



(4) 



where Htq = ea^ + Aa^^, is the Hamiltonian of the trans- 
port qubit (TQ). Here e is the level splitting, and A the 
coherent tunneling within the TQ. H-^^^ = ga^a^a^ is the 
off-resonance dispersive interaction between Dicke sys- 
tem and TQ, Ltq contains the transport properties of 
the TQ^'^, and Lc contains cavity damping terms (e.g., 
photons leaking from the cavity). 
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SLs\^p{t) - 2s\^p{t)sL + p{t)sLs\^ 

s^j^SRp{t) - 2sRp{t)s^j^ + p{t)s^j^SR [5) 
[a^ap — 2apa^ + pa'^a] (6) 



where 



SL = \0){L\, 
SR = \0){R\, 



4 = 



\L){0\, 
\R){0\, 



(7) 
(8) 



Tl and Tr are the left /right tunneling rates for the TQ, 
and 76 is the decay rate of photons out of the cavity 
(throughout, we set = 1). Here p{t) is the density 
matrix describing the state of the qubit-array, cavity, and 
transport qubit system. 



If we assume no back-action from the transport qubit 
onto the Dicke model, the problem is very simple. How- 
ever, the form of the interaction between the transport 
qubit and the effective cavity is still important. As men- 
tioned, off-resonance, e <C CJ, we assume the interac- 
tion is dispersive^, i^int = gcTza^a. For an entirely non- 
destructive passive measurement (with no feedback) , the 
state of the ancillary transport qubit is then just shifted 
by the occupation of the transmission line (i.e. consider- 
ing the mean- field of Eq. [3]) , 



TQ 



(9) 



We are able to calculate the analytical values of {a^ a) 
in the limit N ^ oo. The transport properties are eas- 
ily calculated using a counting-statistics approach, which 
has been well summarised elsewhere^^. Thus, the current 
and zero-frequency current-noise measured through the 
ancillary qubit is simply given by. 
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(11) 



In the limit ^ oo the Dicke Hamiltonian has two 
distinct solutions, corresponding to the two phases of the 
transition. In the superradiant phase both cavity and 
qubit array have a macroscopic mean field displacement. 

In the lower, 'normal phase', we define the occupation 
of the cavity {a^ a) by an effective temperature T and 
frequency Vt^ 
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Where Vl and T depend on the eigenenergies of R'. 



= \ + 'A ± a/K' - ^')' + 16A2u;u;o^l3) 

where e_ is only real for A < Ac, giving the range of this 
solution. The dependence of T and on the eigenvalues 



is via the relations, 
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(14) 
(15) 

(16) 
(17) 
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where P = l/ksT. These define two equations linking 
the three parameters of the cavity /qubit system a;, a;o, 
A, and the three effective parameters of a thermal oscil- 
lator /3, 1], m. By setting one energy scale of the original 
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system such that cj = 1, and that of the thermal oscihator 
such that m = 1, we can uniquely define the correspon- 
dence between the two systems. We use the relations, 



D 

coth(/31]/2) 



l + 2e_e+/L>, 
[sc(e_ -e+)]^ 
L)/(e_s2 + e+c2), 
2e_e+(l + e_e+/L>) 



[cosh{(3n) + l]/sinh(/3cj) 



to obtain, 



c<; me_e+ 



(19) 
(20) 
(21) 

(22) 

(23) 



(24) 



Thus, in this passive measurement regime, in the large 
limit, the occupation of the bosonic mode (which is 
an order parameter of the phase transition) diverges as 
e_ ^ when A ^ Ac- In the next section we discuss the 
effect of this on the current-measurement. 



VI. POWER-LAW SCALING IN TRANSPORT 
PROPERTIES. 

A. Results 

We plot the current and current-noise in Figs. [2] and 
m We immediately see that, at the critical point Ac, 
the large occupation of the cavity mode (which is pro- 
portional to the number of qubits N) acts to blockade 
the current fiow (by "pushing apart" the internal en- 
ergy levels of the transport qubit). Similarly, the zero- 
frequency noise becomes strictly Poissonian at the criti- 
cal point. This is a consequence of the slow current and 
charge-dominated dynamics. Thus, both the current and 
current-noise are operating as signatures, or order pa- 
rameters, of the phase transition^ because of their direct 
dependence on {a^a). 

As mentioned earlier, in previous work^^^ the phase 
transition was studied as a function of multi-qubit- 
oscillator coupling A. However, the transition can also 
be observed for a given constant A, by tuning the en- 
ergy level of the qubits ujo- This is a more realistic ap- 
proach with superconducting qubits as a possible real- 
isation. Qualitatively, the properties of the transition 
are the same. For instance, for A = O.lcc; the transi- 
tion occurs when uJo,c 0.04cj. The sub-radiant phase 
occurs for coq > coq^c^ while the super-radiant phase ap- 
pears when loq < coq^c^ both of which are experimentally- 
accessible regimes. However, because the interaction is 
off-resonance, the convergence to the correct scaling be- 
haviour requires much larger N. 



B. Scaling with the number (N) of qubits 

To observe power-law scaling with we must look 
at the derivative of both the current and current-noise 
with respect to the Dicke multi-qubit-oscillator coupling 
A. The minimum value of these derivatives will act as 
a signature of "precursor behaviour" , and from them we 
can extract the power-law dependence. In Fig. [2jb) we 
show the derivative of the current, and in Fig. [31(a) we see 
that the position of the minimum of the current deriva- 
tive scales as a power law in N via 



(A^ - Ac) (X N- 



-0.68±0.05 



(25) 



This matches a previous result for the scaling of the en- 
tanglement entrop}^^. Similarly the value of the current 
at this minimum point scales logarithmically as 



d{I/e) 
dX 



(X (0.81 ± 0.05) log2 N, 



(26) 



as shown in Fig.[3jb). The value of the current derivative 
obeys similar scaling laws. 

Vidal et al^^ studied the scaling, in A^, at the critical 
point Ac of several properties of the Dicke model. They 
predicted a scaling exponent for l/{a^a)'^ of a = 2/3. 
These exponents are different from those we observe here, 
as they describe behaviour of quantities measured exactly 
at Ac- To extract the same exponents from our numerics 
would require very large values of N. However a recent 
numerical study by Chen et al^^ describes a scheme where 
such exponents can be calculated efficiently for large 
and confirmed^^ the correct exponents for some of these 
quantities. 



VII. BACK-ACTION AND CAVITY LOSS. 

To take into account both the back-action of the trans- 
port qubit, and the loss of photons from the cavity due 
to coupling to the environment, we must solve the entire 
master equation numerically. This is a non-trivial task, 
even with state-of-the-art numerics, and requires careful 
use of sparse-matrix techniques to increase efficiency. 

Dimer et al^^ investigated the thermodynamic limit of 
the Dicke model including losses from the bosonic cavity. 
They found that the critical point was shifted from its 
normal position as a function of the cavity loss 75. In 
Figs. 2(c) and 2(d) we do the same for the finite- A/" case, 
comparing the three possible regimes: zero back-action 
and no cavity loss, zero back-action with cavity loss, and 
a full treatment of cavity loss and back-action. 

In figure 2(c) we see that around the critical point the 
occupancy of the bosonic cavity is almost exactly the 
same for both master equation treatments, but differs 
slightly from the ground state Dicke case. Furthermore 
the strong coupling limit for the full master equation 
treatment saturates because of the bosonic Hilbert space 
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FIG. 2: (Color online) (a) The current I/e versus multi- 
qubit-oscillator coupling A through the transport qubit for 
Tc = 0.1, Tl = Tr = 0.1, e = 0, = cJb = 1, ^ = 0.1 
for iV = 4, 8, 16, 20, 24, oo. (b) The derivative of the current 
through the transport qubit for the same parameter set, ver- 
sus A. Figures (c) and (d) show one particular data curve 
(A/" = 4, 7b = 0.1) for the bosonic occupancy (a^a) and the 
current I/e for the three different approximations; zero back- 
action (ground state of the pure Dicke model), master equa- 
tion with cavity damping, and master equation with cavity 
damping and transport qubit feedback. 
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FIG. 3: (Color online) (a) shows the scaling with N and scal- 
ing exponent of the position (Am) of the minimum of the cur- 
rent derivative: (Am - Ac) oc iy-o-^^^o o^ . Figure (b) shows 
the scaling of the value of the current at this minimum point 
to be (^^^)^ oc (0.81 ± 0.05) logs ^- The parameters used 
here are Tc = 0.1, Fl = F^ = 0.1, e = 0, u ■ 
with data taken at = 4, 8, 16, 20, 24, 40, 60. 



1,^ = 0.1 



cut-oflF needed in solving this complex problem. Further- 
more, in figure 2(d) we see that the full treatment of the 
combined transport qubit/Dicke model shows a reduced 
current profile compared to the two situations with zero 
back-action. This is also the case for other values of N. 

However the coupling to the qubit, and the loss of en- 
ergy from the cavity, has less obvious effects on the prop- 
erties of the phase transition itself. In particular, the 
parity. 



n = exp[z7r(a^a + + j) 



(27) 



is no longer conserved, and the steady state will contain 
components of both the ground state and excited states of 
H]j. Because of this, and the restrictions on the number 
of spins we can efficiently model, it is not possible to 
extract exponents from this data. However, we expect 
the large- limit to still exhibit features of the phase 
transition, as predicted by Dimer et al2^. 



VIII. SIGNATURES OF QUANTUM CHAOS. 

Quantum chaos is a characteristic of non-integrable 
quantum systems. Emary et aP extensively studied the 
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FIG. 4: (Color online) (a) The current-noise F(0) = S(0)/2eI 
versus multi-qubit-oscillator coupling A for: Tc = 0.1, Fl = 
Tr = 0.01, e = 0, = cjb = 1, ^ = 0.1 and for N = 
4,8,16,20,24,00. (b) The derivative versus A. The 

peak scales as a power law of similar to the minimum of 
the current derivative. 
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FIG. 5: (Color online) The inset shows an example of the 
positive imaginary components of the eigenvalues (^f ) of the 
Liouvillian for the damped Dicke model and the main figure 
shows their probability distribution P{E^), normalised to the 
maximum energy gap 5'max, for = 6, A = Ac, and 7?, — 0.1. 
While this contains every possible eigenvalue seperation of the 
Hamiltonian (up to numerical bosonic cut-off), and has not 
been unfolded to remove secular variations, level repulsion is 
still visible. 



(closed) Dicke model and its chaotic properties. In the 
finite- regime they showed that the eigenvalue spec- 
trum of the Dicke model fitted that of the Wigner-Dyson 
distribution^^ when the qubit-boson coupling was around 
the critical point A ^ Ac- Thus, the chaotic behaviour 
is understood to be a 'precursor' of the phase transition, 
driven by the parity conservation at the critical point. 

Here we extend their work by identifying similar dis- 
tributions in the eigenvalues. 



(28) 



of the Liouvillian L which include imaginary components 
i{Ef") from Hd^ as well as real components Ui from the 
cavity loss terms. Here we ignore the back-action and 
electron transport terms here, and focus on the effect of 
cavity damping on the level statistics. 

For the pure-state case (no cavity losses), the von Neu- 
mann equation of motion. 



(29) 



can be written as a set of Nfj coupled equations of the 
matrix elements of p, where Nh is the dimension of the 
Hilbert space for the system described by the Hamilto- 
nian H. If H has Nh eigenvalues E^^ k = 1, A/"//, 



and we take matrix elements according to the eigenba- 
sis of then we can write these linear equations as a 
diagonal matrix with N'jj imaginary eigenvalues 



Et 



--jxk 



Nh 

k,j=i 



(30) 



Every possible energy gap (not just nearest neighbor) in 
the spectrum of H has an eigenvalue in L. 

In Fig. (4) we show the positive branch of the imag- 
inary components of the eigenvalues of L for A/" = 6, 
A = Ac, after removal of the Nh zeros, i.e. the stationary 
states, and the probability distributions of these com- 
ponents. Even though it is not possible to unfold this 
spectrum, and all possible level spacings are present, still 
we see some characteristics of the 'picket-fence' distribu- 
tion^^ of the Rabi Hamiltonian and the universal Wigner- 
Dyson distribution^^. We point out that the eigenvalues 
of this matrix, which is a particular representation of 
the superoperator L, determine many of the higher-order 
transport properties, like the frequency dependant noise. 
This is also seen in scattering theory^^. Further analyt- 
ical work needs to be done to make a strong connection 
between measurable transport quantities, random matrix 
theory, and quantum chaos. 
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IX. FROM CIRCUIT QED TO THE DICKE 
MODEL. 

Al-saidi and Stroud^ have studied a realization of the 
Dicke model using Josephson Junctions coupled to an 
electromagnetic cavity. Operating in the regime 'be- 
tween' charge and flux qubits they showed that, given 
the right parameters, the higher-lying levels of each junc- 
tion can be neglected. In the same way, it is possi- 
ble to derive the Dicke Hamiltonian, Eq. ([1]), from the 
Hamiltonian describing superconducting qubits interact- 
ing with a cavity. The proposal and realization of cavity 
QED^ ^-^-^1-^^1-^^1-^^1-^^1-^^ in a circuit was an important devel- 
opment for quantum optics and condensed matter, and 
thus the observation of strong many-body effects in these 
systems is a natural extension of previous work. 

Alternatively, a large number of qubits, in the form of 
two- level atoms in an atom-chip, coupled to a transmis- 
sion line, was recently proposed as a way to realise the 
large- Dicke model^-. 

X. CONCLUSIONS 

In conclusion, we have shown that current and current- 
noise measurements could be used to test for criticality in 



an 'on chip' experiment. We extracted scaling exponents 
for the Dicke phase transition from semi-analytical and 
numerical modelling, and illustrated how quantum chaos, 
a precursor behaviour to the phase transition, is retained 
in an open-system environment. 
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